arXiv: 1501.04316v2 [hep-ex] 14 Nov 2015 


The transverse momentum dependence of charged kaon Bose-Einstein 
correlations in the SELEX experiment 


The SELEX Collaboration 


G.A. Nigmatkulov 1 , A.K. Ponosov' 1 , U. Akgun", G. Alkhazox3, J. Amaro-Reyes 1 , A. Asratyan 1 , A.G. Atamantchouk 1,1 , A.S. Ayan n , M.Y. Balatzf’ 1 , 
A. Blanco-Covarrubias 1 , N.F. Bondar', P.S. Cooped, L.J. Dauwe 0,1 , G.V. Davidenkcf, U. Dersch 8,3 , A.G. DolgolenkJ, G.B. Dzyubenkcf 1 , 

R. Edelstein b , L. Emediato q , A.M.F. Endler 6 , J. Engelfried 1 , I. Eschrich 8,2 , C.O. Escobar^ 4 , N. Estrada l , A.V. Evdokimovf, I.S. Filimonov 11,1 , 

A. Flores-Castillo 1 , F.G. Garcia q,d , I. Giller k , V.L. Golovtsov 1 , P. Gouffon q , E. GUlmef, M. IorP, S.Y. Jun b , M. Kaya 11,5 , J. Kilmer 4 , V.T. Kim!, 
L.M. Kochenda 1 , I. Konorov 8 ’ 6 , A.P. Kozhevnikov 6 , A.G. Krivshich!, H. Kruger 8,7 , M.A. KubantseJ, V.P. Kubarovsky 6 , A.I. Kulyavtsev 44 , 

N.P. Kuropatkin! ,d , V.F. Kurshetsov 6 , A. Kushnirenko b,e , J. Lach d , L.G. Landsberg 6,1 ,1. Larin 1 , E.M. Leikin b , G. Lopez-Hinojosa!, T. Lungov q , 
V.P. Maleev 1 , D. Mao b,s , P. Mathew 4,9 , M. Mattson b , V. Matveevf, E. McClimenf, M.A. Moines ter 4 , V.V. Molchanov 6 , A. Morelos 1 , A.V. Nemitkin h , 
P.V. NeoustroeV, C. Newsom 11 , A.P. Nilo'J’ 1 , S.B. Nurushev 6 , A. Ocherashvili 4,10 , Y. Onel n , S. Ozkorucuklu' 1,11 , A. Penzo r , S.V. Petrenko 6 , 

M. Procario 4,12 , V.A. Prutskoi 1 , B.V. Razmyslovich/’ 13 , D.A. Romanov 1 , V.I. Rud h , J. Russ 4 , J.L. Sanchez-Lopez 1 , A.A. Savchenko 1 , J. Simon 8,14 , 
G.V. Sinev 1,15 , A.I. Sitnikox$, V.J. Smith! 11 , M. Srivastava q , V. Steiner k , V. Stepanov 1,13 , L. Stutte d , M. Svoiski 1,13 , V.V. Tarasovf, N.K Terentyev 1,4 , 

/. Torres 1,16 , L.N. Uvarov 1 , A.N. Vasiliev 6 , D.V. Vavilov 6 , E. Vazquez-Jduregui 1,17 , V.S. Verebryusov 4,1 , V.A. Victorov 6 , V.E. Vishnyakovf, 

A.A. Vorobyov 1 , K. Vorwalter 8,18 , J. You b,d , R. Zukanovich-Funchal q 

a Bogazici University, Bebek 80815 Istanbul, Turkey 
4 Carnegie-Mellon University, Pittsburgh, PA 15213, U.S.A. 

6 Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, Brazil 
d Fermi National Accelerator Laboratory, Batavia, IL 60510, U.S.A. 

6 Institute for High Energy Physics, Protvino, Russia 
* Institute of Theoretical and Experimental Physics, Moscow, Russia 
8 Max-Planck-Institut fiir Kernphysik, 69117 Heidelberg, Germany 
4 Moscow State University, Moscow, Russia 

1 National Research Nuclear University MEPhI (Moscow Engineering Physics Institute), Moscow, Russia 
1 Petersburg Nuclear Physics Institute, St. Petersburg, Russia 
k TelAviv University, 69978 Ramat Aviv, Israel 
1 Universidad Autonoma de San Luis Potosi, San Luis Potosi, Mexico 
m University of Bristol, Bristol BS8 1TL, United Kingdom 
n University of Iowa, Iowa City, 1A 52242, U.S.A. 

° University of Michigan-Flint, Flint, Ml 48502, U.S.A. 
p University of Rome “La Sapienza” and INFN, Rome, Italy 
q University ofSao Paulo, Sao Paulo, Brazil 
r University of Trieste and INFN, Trieste, Italy 


Abstract 

We report on the measurement of the one-dimensional charged kaon correlation functions using 600 GeV/c 2T, n 
and 540 GeV/c p beams from the SELEX (E781) experiment at the Fermilab Tevatron. K L correlation functions 
are studied for three transverse pair momentum, kj, ranges and parameterized by a Gaussian form. The emission 
source radii, R, and the correlation strength. A, are extracted. The analysis shows a decrease of the source radii with 
increasing kaon transverse pair momentum for all beam types. 
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1. Introduction 


In this paper we present results from K L K L correla¬ 
tion femtoscopy study in 600 GeV/c YrC(Cu), n C(Cu) 
and 540 GeV/c pC{Cu) interactions from the SE- 
LEX (E781) experiment [j7|] at the Fermilab Tevatron. 
The correlation femtoscopy method allows to study 
spatio-temporal characteristics of the emission source 
at the level of 1 fm = 10~ 15 m. The method is based 
on the Bose-Einstein enhancement of identical boson 
production at small relative momentum. The quantum 
statistics correlations were first observed as an enhanced 
production of the identical pion pairs with small open¬ 
ing angles in proton-antiproton collisions 0. In 1960 
this enhancement was explained by the symmetriza- 
tion of the two-particle wave function by G. Goldhaber, 
S. Goldhaber, W.-Y. Lee, and A. Pais (GGLP effect) 0]. 
Later, in the 1970s, Kopylov and Podgoretsky suggested 
studying the interference effect in terms of the correla¬ 
tion function. They proposed the mixing technique to 
construct the uncorrelated reference sample, and clari¬ 
fied the role of the space-time characteristics of parti¬ 
cle production Subsequently, two-particle cor¬ 

relations at small relative momentum were systemati¬ 
cally studied for lepton-lepton 0, lepton-hadron 0, 
hadron-hadron 0, and heavy-ion HpUlI l collisions. It 
was found that the system created in heavy-ion colli¬ 
sions undergoes the collective expansion and may be de¬ 
scribed by relativistic fluid dynamics 1112- 11511 . By using 
the width of the quantum statistical enhancement, one 
can measure the radii R of the emitting source. The de¬ 
crease of the extracted radii with increasing transverse 
pair momentum may be interpreted as the decrease of 
the “homogeneity lengths” 1161] due to collective trans¬ 
verse flow. 

A comparison of femtoscopic measurements in lep¬ 
ton and hadron-induced 117111811 collisions with heavy- 
ion collisions shows similar systematics [19, 20]. These 
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studies usually performed for pions. However, mea¬ 
surements of heavier particles may provide additional 
information about the size, orientation and dynamical 
timescales of the emission region. 


2. Experimental setup and data selection 


The SELEX (E781) detector is a three-stage magnetic 
spectrometer designed for charm hadroproduction study 
at xf >0.1 (xp = ^—). We report on the analysis of 
1 billion events of 600 GeV/c E~C(Cm), n~C{Cu ) and 
540 GeV/c p C(Cu) interactions recorded during the 
1996-1997 fixed target run. About 2/3, 1/6 and 1/6 of 
the data were obtained on IT, n and p beams, respec¬ 
tively. 

The beam particle was identified as a meson or 
a baryon by a transition radiation detector. Interac¬ 
tions occurred on segmented targets, which consisted 
of 2 copper and 3 diamond foils separated by 1.5 cm 
clearance, and had a total thickness of 5% of an inter¬ 
action length for protons. Particles were tracked in a 
set of 20 vertex Silicon Strip Detectors (SSD) arranged 
in 4 sets of planes with a strip pitch of 20-25 pm, ro¬ 
tated by 45°. Each of the detectors has a hit detection 
efficiency greater than 98%. Transverse vertex position 
resolution (cr) was 4 pm for the 600 GeV/c beam tracks. 
The average longitudinal vertex position resolution was 
270 pm for primary vertex and 560 pm for secondary 
vertex. The detector geometry covers the forward 150 
mrad cone. The particle momentum was measured by 
deflection of the track position by two magnets Ml and 
M2 in a system of proportional wire chambers and sil¬ 
icon strip detectors. Momentum resolution of a typical 
100 GeV/c track was cr p lp ~ 0.5%. A Ring Imaging 
Cherenkov detector (RICH) performed particle identi¬ 
fication in a wide momentum range and provided 2 cr 
Kjn separation up to 165 GeV/c and single track ring 
radius resolution of 1.4% |2lj]. The kaon identifica¬ 
tion efficiency was over 90% above the kaon threshold 
(«43 GeV/c). The average number of tracks reaching 
the RICH was about 5 per event IL22ll. The layout of the 
spectrometer is described elsewhere 0. 

In this analysis we used primary tracks that have 
vertex silicon track segment matched with downstream 
segments measured in the Ml and M2 spectrometers, 
with the momentum from 45 to 160 GeV/c. In order to 
reduce the contamination of secondary particles, it was 
required that the extrapolated track distance to the pri¬ 
mary vertex was less than 15 pm in the transverse plane. 
Only tracks that matched the RICH detector were used 
in the analysis. Charged kaons were identified with the 
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likelihood to be a kaon at least three times exceeding 
any other particle hypothesis. Fig. |TJa) shows the sin¬ 
gle kaon purity as a function of the momentum for the 
E n and p beams. It is defined as the fraction of the 
accepted kaon tracks that correspond to true kaon parti¬ 
cles. The single particle purity was estimated from the 
RICH ring radius distributions of the data and by study¬ 
ing PYTHIA [23h simulations with the particle embed¬ 
ding through the detector. The main contamination for 
charged kaons in the momentum range p > 120 GeV/c 
comes from pions, because of the close ring radii in the 
RICH detector. 
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average momentum K = {p\ + /A)/2 of two particles: 


c(4J) = ^^-d(4,£), 

B(q, K) 


( 2 ) 


where A(q, K) is the measured distribution of relative 
momentum within the same event, B(q, K) is the ref¬ 
erence or background distribution that is similar to the 
experimental distribution in all respects except for the 
presence of femtoscopic correlations. The reference 
sample is usually formed from particles that come from 
different events (event mixing technique Jbtl j. The quan¬ 
tity D(q, K) is a so-called correlation baseline that de¬ 
scribes all non-femtoscopic correlations, such as, for 
instance, the correlations caused by the energy and mo¬ 
mentum conservation-induced correlations [[24]. In or¬ 
der to eliminate possible biases due to the construction 
of the reference samples, the measured correlation func¬ 
tions were corrected on the simulated distributions by 
constructing the double ratio: 


O - (dN k+ k+ Id.Q\ IdN mc^k*I dQ\ 
J ~ ! dN„f/dQ r\ dN MC ,refldQ )’ 


(3) 


Figure 1: Single S'* purity as a function of momentum (a) and the 
average transverse pair momentum dependence of the S'* pair purity 
for the X (circles), it (squares) and p (crosses) beams (b). 


The electrons were eliminated from the analysis us¬ 
ing the transition radiation detector (ETRD) that was 
placed before RICH. The contamination from other par¬ 
ticle species in the studied momentum range is negligi¬ 
ble. Fig. |TJb) shows the charged kaon pair purity as a 
function of the average transverse momentum of the pair 
obtained for the E , n~ and p beams. The K L pair purity 
is calculated as a product of two single-particle purities 
using the experimental momentum distributions. 

After applying the cuts 4842147,597101 and 103551 
identical charged kaon pairs were selected for E~, n~ 
and p beams, respectively. 


3. Correlation femtoscopy 


The two-particle correlation function is defined as the 
ratio of the probability to measure two particles with 
momenta p\ and /A to their single particle probabilities: 


C(pi,pi) 


P(Pl,P2) 
P(p\)P(pi) 


( 1 ) 


Experimentally, one studies the correlation function 
C(q, K) in terms of relative momentum q — p\ - /A and 


where the subscripts “MC” and “MC, re/” correspond 
to the simulated data. 

By virtue of the limited statistics available for 
the 7r“ and p beams, only the one-dimensional fem¬ 
toscopic charged kaon analysis of correlation func¬ 
tions in terms of invariant relative momentum, Q = 
V(/?i - /A) 2 - ( E\ - E 2 ) 2 , was performed. In order to 
extract the size of the emission region, R, one can use 
the Goldhaber parametrization. This assumes that the 
emitting source of identical bosons is described by a 
spherical Gaussian density function: 


C(Q) = N(l+Ae- R2Q2 )-D(Q), 


(4) 


where A is a normalization factor, A describes the cor¬ 
relation strength, and D(Q) is the baseline distribution. 
In the current analysis, the second order polynomial, 
D(Q) - 1 + aQ + bQ 2 , was used for estimation of the 
baseline distribution. The momentum correlations of 
particles emitted at nuclear distances are also influenced 
by the effect of final-state interaction (FSI), Coulomb 
and strong interactions 125 - 28|]. For identical kaons, the 


effect of strong interactions is negligible [29]. The cor¬ 
relation function of identical bosons should increase at 
low relative momentum, except for small values where 
Coulomb interaction becomes dominant. This may be 
taken into account by modifying Eq. ©: 


C(Q) = N ((1 - T) + AK{Q) (l + e- RlQ2 ))D(Q), (5) 
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where the factor K(Q) is the squared like-sign kaon 
pair Coulomb wave function integrated over a spherical 
Gaussian source SQ. 


4. Results and discussions 

The results discussed in this Letter were obtained 
with the same detector setup, cuts and fitting proce¬ 
dures, giving an opportunity to compare the properties 
of the emission region for different hadron-induced col¬ 
lisions. The analysis was performed for three average 
transverse pair momentum kp - \pri + Pti\ /2 ranges: 
(0.00 - 0.30), (0.30 - 0.55), (0.55 - 1.00) GeV/c and for 
the three beam types: X n , p. The event mixing tech¬ 
nique was used to construct the uncorrelated reference 
sample. Only events with two or more identical charged 
kaons, grouped by production target, were used in the 
event mixing. Kaons from adjacent events for each tar¬ 
get were combined to provide an uncorrelated exper¬ 
imental background. Due to small differences in the 
measured correlation functions, the positive and nega¬ 
tive kaon four-momentum distributions were combined 
in the numerator and the denominator before construct¬ 
ing the ratio. A purity correction was applied to the 
experimental correlation functions according to the ex¬ 
pression: 

Cexperimentali.Q') ~ 1 


C{Q) 


P(Q ) 


+ 1 , 


( 6 ) 


where P(Q ) is the pair purity. 

The top row of the Fig. [2] shows the experimental cor¬ 
relation functions (solid circles) after applying the pu¬ 
rity correction measured in the 0.30 < kj < 0.55 GeV/c 
region for IT, n~ and p beams. The correlation func¬ 
tions were normalized such that C(Q ) = 1 for 0.5 < 
Q < 0.7 GeV/c, where Bose-Einstein correlations are 
absent and the influence of the non-femtoscopic effects 
is small. The deviation of the correlation functions 
from unity at Q > 0.7 GeV/c correspond to the non- 
femtoscopic correlations and the well-defined enhance¬ 
ment at Q < 0.4 GeV/c is due to the quantum statistical 
correlations. The Coulomb repulsion between like-sign 
kaons leads to the decrease of the correlation functions 
atg < 0.1 GeV/c. 

In order to correct for non-femtoscopic effects, the 
Monte Carlo event generator PYTHIA-6.4.28 [23] with 
different tunes (Perugia 0, Perugia 2010 and Peru¬ 
gia 2011 l33ll ) was used. PYTHIA contains neither 
Bose-Einstein correlations nor the final-state interac¬ 
tions. On the other hand, PYTHIA contains other 
kinematic effects, for instance, energy and momen¬ 
tum conservation effects, that could lead to baseline 


correlations. The Perugia 2011 tune, which best de¬ 
scribes charged-particle multiplicity, was used as the 
main minimum-bias MC sample. The top row of the 
Fig. El shows the comparison of simulated correlation 
functions (empty circles), where PYTHIA events were 
filtered through the analysis cuts, with the experimental 
distributions (solid circles) measured for 0.30 < kj < 
0.55 GeV/c range. The PYTHIA-generated correlation 
functions were normalized in the same way as the ex¬ 
perimental correlation functions. 
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Figure 2: Top row shows the experimental (full circles) and PYTHIA- 
generated (open circles) correlation functions of identical kaons ob¬ 
tained from using the event mixed reference samples. The bottom 
row shows the fits to the double ratios according to Eq. {3). The sec¬ 
ond order polynomial was used for estimation of the non-femtoscopic 
effects. The correlation functions are measured for 0.30 < kp < 
0.55 GeV/c range and columns from left to right represent the data 
obtained for X . jt and p beams, respectively. 


It is seen that PYTHIA qualitatively describes the 
experimental baseline in the region Q > 0.5 GeV/c, 
where the effect of femtoscopic correlations is negligi¬ 
ble. Since the MC calculation does not include wave 
function symmetrization for identical particles, the fem¬ 
toscopic peak at low relative four-momentum region, 
Q < 0.4 GeV/c, is absent. 

The bottom row of the Fig. [2] shows double ratios, 
where the experimental correlation functions are di¬ 
vided by the PYTHIA-generated ones, obtained for the 
X - , 7r - and p beams in the 0.30 < k r < 0.55 GeV/c re¬ 
gion. The double ratios were fitted using Eq. ([5]). In the 
current analysis, the Coulomb function K(Q) was inte¬ 
grated over a spherical source of 1 fm. Due to imperfec¬ 
tions of the simulation in the Q > 0.7 GeV/c region, the 
non-femtoscopic term D(Q) — 1 + aQ + bQ 2 was used. 

To estimate the influence of choice of the reference 
sample, the different methods of constructing uncorre¬ 
lated charged particle distributions were used: opposite- 
charge pairs and rotated particles, where pairs are con¬ 
structed after inverting the x and y components of the 
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Table 1: Systematic uncertainty (minimal and maximal) values for different sources of systematic uncertainty (in percent). The values of minimum 
(maximum) uncertainty can be from different average transverse pair momentum range or the beam type, but from a specific source. 


The systematic uncertainty source 

A(%) 

R(%) 

Event/particle selection 

1-7 

1-9 

PID and purity 

0-4 

0-6 

Fit range 

1-5 

1-4 

Momentum resolution 

0-1 

0-1 

Two-track effects 

- 

- 

Non-femtoscopic form 

0-4 

1-11 

Coulomb function 

- 

- 

Reference sample 

1-8 

5-13 

Total 

2-13 

5-21 


three-momentum of one of the two particles. Fig. [3] 
shows the double ratios obtained from using rotated 
particles (top panel) and opposite-charge pairs (bottom 
panel) reference samples. The double ratios were fit¬ 
ted using Eq. (0; and the second order polynomial, 
D(Q) - 1 + aQ + bQ 2 , was used to describe the non- 
femtoscopic term. It was found that the extracted femto- 
scopic parameters, A and R, obtained from using rotated 
particles reference samples are similar to those from the 
event mixed ones. 
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event mixed ones. The range 0.19 < Q < 0.35 GeV/c 
was excluded from fits due to the 0(1020) meson decay 
and because the influence of the final-state interactions 
between opposite-charge kaon pairs was not taken into 
account. 

The different sources of systematic uncertainties were 
studied for each k T range and beam type. Tab. [T] shows 
the maximal and minimal values of systematic uncer¬ 
tainty that correspond to specific uncertainty sources. 
The values of the total uncertainty are not necessarily 
equal to the quadratic sum of all the uncertainties due to 
the fact that they can come from different beam types or 
average transverse pair momentum ranges. 

The uncertainty due to the fit range was estimated by 
varying the upper limit of the fit from Q - 1 GeV/c to 
Q = 0.6 GeV/c. The lowest value of the upper limit of 
the fit range corresponds to the end of the correlation 
region. Changing the radius of the Coulomb source in 
the range from 0.5 fm to 1.5 fm has negligible effect on 
the extracted emitting source parameters. 

Different baseline shapes 134 - 3f| were used to esti¬ 
mate the systematic uncertainty due to the baseline de¬ 
termination: 


Figure 3: The correlation functions constructed with rotated particles 
(top panel) and opposite-charge pairs (bottom panel) for 0.30 < k? < 
0.55 GeV/c range. The range 0.19 < Q < 0.35 GeV/c on the bottom 
panel is excluded from fits due to the contribution from the 0(1020) 
meson decay. The fits were performed using Eq. j5j. The columns 
represent X (left), n~ (middle) and p (right) beams. 


D(Q) = 1, 

(7) 

0(0 = 1 + aQ, 

(8) 

D(Q)= l+e- aQ \ 

(9) 


The reference samples constructed from the opposite- 
charge kaon pairs contain peaks coming from strong 
resonance decays; and are influenced by the Coulomb 
attraction at Q < 0.1 GeV/c. The magnitudes of the res¬ 
onance peaks measured for the real and simulated corre¬ 
lation functions were found to be different. This can be 
explained by the absence of the final-state rescattering 
of particles in PYTHIA. The double ratios obtained for 
unlike-sign kaon pairs were fitted the same way as the 


D(Q) = Vl +aQ 2 + bQ 4 . (10) 

The smearing of single particle momenta was studied 
by embedding simulated kaon tracks with known mo¬ 
menta through the detector. Experimental correlation 
functions were corrected for momentum resolution us¬ 
ing the expression: 

(rx\ Cuncorr ( Q)C unsmeare m .... 

L'COrAtJ) — ^ s sy. ? 

^ smeared ) 
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where C corr (Q ) is the corrected correlation function and 
CuncorriQ ) represents the measured correlation function. 
The C lmsmeared (Q ) and C smeared (Q ) are the correlation 
functions without and with taking into account the mo¬ 
mentum resolution effect, respectively. The smearing 
of the single particle momenta leads to the smearing of 
the correlation function in the Q < 0.2 GeV/c range. 
This decreases the amplitude of the peak and makes it 
wider. The systematic uncertainty of this effect does 
not exceed 1%. The two-track reconstruction effects: 
“merging”, when two tracks are reconstructed as one; 
and “splitting”, when one track is reconstructed as two, 
were studied and found to be negligible. 

The main systematic uncertainty arises from using 
different methods of constructing the uncorrelated ref¬ 
erence sample. The dominant contribution comes from 
using opposite-charge kaon pairs. The contamination of 
the “Reference sample” uncertainty from using different 
PYTHIA tunes does not exceed 5%. 

The systematic errors on R and A for each beam type 
and kj range are taken as the rms spread of the values 
obtained for the different sources of systematic uncer¬ 
tainty. 

The results of fits of R(Q ) based on the parametriza- 
tion of Eq. (0 with D(Q) = 1 + aQ + bQ 2 are given 
in Table [2] The extracted source radii and correla¬ 
tion strength for X (circles), n~ (squares) and p (stars) 
beams as a function of transverse pair momentum are 
shown in Figs. [4]' a) and [4] b), respectively. The source 
radii slightly decrease with increasing k T for all the 
beam types, except the highest average transverse pair 
momentum interval for the n beam. The femtoscopic 
radii measured for X , n~ and p beams are consistent 
within the uncertainties. The small difference between 
measured source parameters probably arises from dif¬ 
ferent contamination from resonance decays & 

The decrease of the source radii with transverse pair 
momentum was previously observed in heavy-ion colli¬ 
sions and interpreted as a collective hydrodynamic be¬ 
havior (collective flow) |3f| 39]. The first direct com¬ 
parison of correlation femtoscopy in p + p and heavy- 
ion collisions under the same detector conditions, re¬ 
construction, analysis and fitting procedures was per¬ 
formed by the STAR collaboration |2Ctl . It was shown 
that p + p collisions also have the transverse momentum 
scaling. Although the interpretation of these results is 
still unclear, the similarities could indicate a connection 
between the underlying physics. 

The transverse momentum dependence was also ob¬ 
served for nn in e + e [HqI [ill] and pp collisions ll20[ l35. 


4211 . For the first time a similar analysis of charged kaon 


Bose-Einstein correlations for more than one transverse 



Figure 4: source parameters R (a) and A (b) measured for X 

(circles), n (squares) and p (stars) beams as a function of transverse 
pair momentum, k j . Statistical (lines) and systematic (boxes) uncer¬ 
tainties are shown. 


pair momentum and multiplicity range was recently per¬ 
formed by the ALICE collaboration at the LHC in pp 
collisions at yfs = 7 TeV [34]. It was shown that 
charged kaon femtoscopic radii decrease with trans¬ 
verse pair momentum for middle and high multiplicity 
ranges. 

There are several possible processes that may lead to 
the kj dependencies in the hadronic collisions: 

1. The space-momentum dependence of the fem¬ 
toscopic radii may be generated by long-lived reso¬ 
nances [43]. In particular this may play a significant 
role in high multiplicity bins, where the bulk collective 
flow is predicted [44]. 

2. Humanic’s model [45]], based on space-time ge¬ 
ometry of hadronization and effects of final-state rescat¬ 
tering between hadrons, reproduces both multiplicity 
and transverse mass dependence measured at the Teva- 
tron 1461. 

3. In small systems, the string fragmentation should 
generate momentum and space correlations, such as kj 
dependence of the source radii. However, there are al¬ 
most no quantitative predictions that may be directly 
compared with data except the r-model in which space- 
time and momentum space are strongly correlated 1 4711 . 
Moreover, the Lund string model is not able to repro¬ 
duce the mass dependence of the radii [17, 48 - 50h. 

4. Hydrodynamic bulk collective flow may lead to a 
kj dependence that is very similar to that from heavy- 
ion collisions. 

Taking the aforementioned possibilities, the origin of 
the transverse pair momentum dependence of the fem¬ 
toscopic radii in hadronic collisions is still unclear. Fur¬ 
ther theoretical studies are needed in order to understand 
the underlying physics. 
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Table 2: K' K + source parameters for X , n and p beams. Statistical and systematic uncertainties are presented. 


Beam type 

k T (GeV/c) 

*7AW 

N 

A 

R (fm) 

a (GeV/cT 1 ) 

b (GeV/c’-) 


(0.00 - 0.30) 

126/45 

1.23 ± 0.01 

0.71 ±0.02 ±0.08 

1.32 ±0.02 ±0.07 

-0.59 ± 0.02 

0.38 ± 0.02 

IT 

(0.30-0.55) 

85/45 

1.18 ±0.01 

0.66 ± 0.02 ± 0.08 

1.18 ±0.02 ±0.05 

-0.47 ± 0.03 

0.28 ± 0.02 


(0.55 - 1.00) 

142/45 

1.05 ±0.03 

0.66 ±0.04 ±0.10 

0.98 ± 0.03 ± 0.04 

-0.22 ± 0.08 

0.13 ±0.05 


(0.00 - 0.30) 

62/45 

I.I9 ± 0.03 

0.67 ± 0.06 ± 0.09 

1.25 ±0.06 ±0.06 

-0.49 ± 0.07 

0.31 ±0.05 

n~ 

(0.30-0.55) 

66/45 

1.21 ± 0.04 

0.69 ± 0.06 ± 0.06 

1.13 ±0.06 ±0.06 

-0.46 ± 0.09 

0.24 ± 0.07 


(0.55 - 1.00) 

58/45 

1.34 ± 0.07 

0.44 ±0.10 ±0.11 

1.16 ±0.19 ±0.14 

-0.71 ± 0.09 

0.42 ± 0.09 


(0.00 - 0.30) 

65/45 

1.51 ±0.06 

0.98 ±0.17 ±0.13 

1.54 ±0.16 ±0.17 

-0.97 ±0.10 

0.62 ± 0.08 

p 

(0.30 - 0.55) 

62/45 

1.39 ±0.12 

0.80 ±0.15 ±0.13 

1.32 ±0.12 ±0.15 

-0.72 ±0.13 

0.40 ±0.11 


(0.55 - 1.00) 

43/44 

1.26 ± 0.16 

0.91 ±0.24 ±0.11 

1.13 ±0.17 ±0.11 

-0.61 ±0.31 

0.37 ± 0.24 


5. Summary 

Charged kaon Bose-Einstein correlations were mea¬ 
sured in the SELEX experiment. One-dimensional 
charged kaon correlation functions in terms of the in¬ 
variant four-momentum difference were constructed for 
IT, n~ and p beams and three transverse pair momen¬ 
tum ranges: (0.00 - 0.30), (0.30 - 0.55) and (0.55 - 
1.00) GeV/c. The source parameters of correlation 
strength. A, and source radii, R, were extracted for all 
beam types and for the three average transverse pair mo¬ 
mentum ranges. The slight decrease of the femtoscopic 
radii with pair transverse momentum was observed for 
all three beam types, except for the highest kj range of 
the tC beam. The values of the emitting source radii ob¬ 
tained for S“, 7 r“ and p beams are consistent within the 
uncertainties. 
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